Fragmentation of shells 
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We present a theoretical and experimental study of the fragmentation of closed thin shells made 
of a disordered brittle material. Experiments were performed on brown and white hen egg-shells 
under two different loading conditions: impact with a hard wall and explosion by a combustible 
mixture both give rise to power law fragment size distributions. A three-dimensional discrete element 
model of shells is worked out. Based on simulations of the model we give evidence that power law 
fragment mass distributions arise due to an underlying phase transition which proved to be abrupt 
for explosion and continuous for impact. We demonstrate that the fragmentation of closed shells 
defines a new universality class of fragmentation phenomena. 

PACS numbers: 46.50.-|-a, 62.20.Mk, 64.60.-i 



Closed shells made of solid materials are often used in 
every day life and industrial applications in form of con- 
tainers, pressure vessels or combustion chambers. From 
a structural point of view aircraft vehicles, launch vehi- 
cles like rockets and building blocks of a space station 
are also shell-like systems and even certain types of mod- 
ern buildings can be considered as shells. The egg-shell 
as nature's oldest container proved to be a reliable con- 
struction for protecting life. In many of the applications 
shell-like constructions operate under an internal pres- 
sure and they usually fail due to an excess internal load 
which can arise either as a result of slowly driving the 
system above its stability limit, or by a pressure pulse 
caused by an explosive shock inside the shell. Due to 
the widespread applications, the failure of shell systems 
is a very important scientific and technological problem 
which can also have an enormous social impact due to 
the human costs arising, for instance, in accidents. 

Fragmentation, i.e. the breaking of particulate ma- 
terials into smaller pieces is abundant in nature occur- 
ring on a broad range of length scales from meteor im- 
pacts through geological phenomena and industrial ap- 
plications down to the break-up of large molecules and 
heavy nuclei 0, 0> S I3- The most striking obser- 
vation concerning fragmentation is that the distribution 
of fragment sizes shows a power law behavior indepen- 
dent on the way of imparting energy, relevant micro- 
scopic interactions and length scales involved, with an 
exponent depending only on the dimensionality of the 
system |1I1I1I1I1I1|3,|1I3- Detailed experimental 
and theoretical studies revealed that universality prevails 
for large enough input energies when the system falls 
apart into small enough pieces P,I3,SSI1|> however, at 
lower energies a systematic dependence of the exponent 
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FIG. 1: Time evolution of the explosion of an egg-shell. 



on the input energy is evident P|. Recent investigations 
on the low energy limit of fragmentation suggest that 
the power law distribution of fragment sizes arises due 
to an underlying critical point Former stud- 

ies on fragmentation have been focused on bulk systems 
in one, two and three dimensions, however, hardly any 
studies have been devoted to the fragmentation of shells. 
The peculiarity of the fragmentation of closed shells origi- 
nates from the fact that their local structure is inherently 
two-dimensional, however, the dynamics of the systems, 
the motion of material elements, deformation and stress 
states are three-dimensional which allows for a rich vari- 
ety of failure modes. 

In this letter we present a theoretical and experimental 
study of the fragmentation of closed brittle shells arising 
due to an excess load inside the shell. We performed 
experiments on the explosion and impact fragmentation 
of hen egg-shells resulting in power law fragment size 
distributions. Based on large scale molecular dynam- 
ics simulations of a discrete element model of shells we 
give evidence that power law fragment mass distributions 
arise due to an underlying phase transition which proved 
to be abrupt for explosion and continuous for impact. 
It is shown that the fragmentation of closed shells be- 
longs to a universality class different from that of the 
two- and three-dimensional bulk systems. Our results 
have also important implications for the fragmentation 
of bulk systems. 



For simplicity, our experiments were performed on or- 
dinary brown and white egg-shells which consists of a 
natural disordered material with a high degree of brittle- 
ness. In the preparation, first two holes of regular circu- 
lar shape were drilled in the bottom and top of the egg 
through which the content of the egg was blown-out. The 
inside was carefully washed and rinsed out several times 
and finally the empty shells were dried in a microwave 
oven to get rid of all moisture of the egg-shell. 

In the impact experiments, intact egg shells are cata- 
pulted against the ground at a high speed using a simple 
setup of rubber bands. In explosion experiments, initially 
the egg shell is flooded with hydrogen and hung vertically 
inside a plastic bag. The combustion reaction is initiated 
by igniting the escaping hydrogen on the top of the egg. 
Due to the escaping hydrogen, oxygen is drawn-up into 
the egg through the bottom hole. When enough air has 
entered to form a combustible mixture inside the egg, 
the flame penetrates through the top hole and the egg 
explodes. Both kinds of experiments were carried out in- 
side a soft plastic bag so that secondary fragmentation 
due to the impact of fragments did not occur. 

The resulting egg-shell pieces are then carefully col- 
lected and placed on the tray of a scanner without over- 
lap. In the scanned image fragments appear as black 
spots on a white background, which were further ana- 
lyzed by a cluster searching code. The mass of fragments 
was determined as the number of pixels of spots in the 
scanned image. 

Three consecutive snapshots of the explosion process 
of an egg-shell are presented in Fig. ^ taken by a high 
speed camera of 400 Hz. Based on the snapshots the 
total duration of an explosion is estimated to be of the 
order of 1 msec. An example of scanned egg-shell pieces 
can be seen in the inset of Fig. |5| Since the pressure 
which builds up during combustion can be slightly con- 
trolled by the hole size, i.e. the smaller the hole is, the 
higher the pressure at the explosion is, we performed sev- 
eral series of explosion experiments with hole diameters 
d between 1.2 and 2.5 millimeters. Fig. |21 presents the 
fragment mass distributions F{m) for impact and explo- 
sion experiments averaged over 10-20 egg-shells for each 
curve. For the impact experiment, a power law behav- 
ior of the fragment mass distribution F(rn) ^ can 
be observed over three orders of magnitude where the 
value of the exponent can be determined with high pre- 
cision to T = 1.35 ± 0.02. Explosion experiments result 
also in a power law distribution of the same value of r 
for small fragments with a relatively broad cut-off for 
the large ones. Smaller hole diameter d in Fig. |21 i.e. 
higher pressure, gives rise to a larger number of fragments 
with a smaller cut-off mass and a faster decay of the dis- 
tribution F{m) at the large fragments. Comparing the 
number of fragments obtained, the ratio of the pressure 
values in the explosions at hole diameters d = 1.2 and 
2.0 mm, presented in Fig. |21 was estimated to be about 
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FIG. 2: (Color online) Comparison of fragment mass distri- 
butions obtained in explosion experiments with two hole sizes, 
and in the impact experiment to the simulation results. Inset: 
scanned pieces of an impact experiment. 



1.6. Note that the relatively small value of the exponent 
T can indicate a cleavage mechanism of shell fragmenta- 
tion and is significantly different from the experimental 
and theoretical results on fragmenting two-dimensional 
bulk systems where 1.5 < t < 2 has been found, and 
from the three-dimensional ones where r > 2 is obtained 

Qiiiiiiii. 

Our theoretical study is restricted for simplicity to 
spherical shells, i.e. we worked out a three-dimensional 
discrete element model of spherical shells by discretiz- 
ing the surface of the unit sphere into randomly shaped 
triangles (Delaunay triangulation) . The nodes of the tri- 
angulation represent point-like material elements in the 
model whose masses are defined by the area of the dual 
Voronoi polygon assigned to it. The bonds between 
nodes are assumed to be springs having linear elastic be- 
havior up to failure. Disorder is introduced in the model 
solely by the randomness of the tessellation so that the 
masses of the nodes, the lengths and the cross-sections of 
the springs are determined by the tessellation (quenched 
structural disorder). After prescribing the initial con- 
ditions of a specific fragmentation process studied, the 
time evolution of the system is followed by solving the 
equation of motion of the nodes. In order to account 
for crack formation in the model, springs are assumed 
to break during the time evolution of the system when 
their deformation e exceeds a fixed breaking threshold Sc 
resulting in a random sequence of breakings due to the 
disordered spring properties. As a result of successive 
spring breakings cracks nucleate, grow and merge on the 
spherical surface giving rise to the complete break-up of 
the shell into pieces. Fragments of the shell are defined 
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FIG. 3: Final fragmented states for a) impact and b) pressure 
pulse loading. Particle positions are projected back to the 
initial state on the surface. Black lines indicate cracks. 



in the model as sets of nodes connected by the remaining 
intact springs. The process is stopped when the system 
has attained a relaxed state. 

In the computer simulations two different ways of load- 
ing have been considered which model the experimental 
conditions and represent limiting cases of energy input 
rates: (i) pressure pulse and (ii) impact load starting from 
an initially stress free state. A pressure pulse in a shell 
is carried out by imposing a fixed internal pressure Pq 
from which the forces acting on the material elements 
are calculated. The constant pressure Pq gives rise to 
an expansion of the system with a continuous increase of 
the imparted energy Etot = Po^V, where AV denotes 
the volume change with respect to the initial volume Vo- 
Since the force F acting on the shell is proportional to 
the actual surface area A, the system is driven by an in- 
creasing force F ^ Po ■ A during the expansion process. 
The impact loading realizes the limiting case of instan- 
taneous energy input Etot = Eg by giving a fixed ini- 
tial radially oriented velocity Vo to all material elements. 
Simulations performed varying the control parameters Pq 
and Eg in a broad range revealed that for both loading 
cases there exists a critical value Pc and Ec below which 
the shell keeps its integrity suffering only partial failure 
in the form of cracks (damaged state), above it however, 
complete break-up occurs into pieces (fragmented state) . 
Since the break-up of the shell under pressure (impact) 
loading is analogous to the stress (strain) controlled frac- 
ture of disordered bulk solids, the transitions between the 
damaged and fragmented states arises abruptly (contin- 
uously) at the critical point. 

Quantitative characterization of the break-up process 
when increasing the control parameter can be given by 
monitoring the average mass of the largest fragment nor- 
mahzed by the total mass {M^ax / Mtot) ■ M^ax is a 
monotonically decreasing function of both Pg and Eq, 
however, the functional forms are different in the two 
cases, see Fig.^ At low pressure values in Fig.^M„iax 
is practically equal to the total mass since hardly any 
fragments are formed. Above Pc however, M^ax gets sig- 
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FIG. 4: {Mmax/Mtot) versus Po and Eo. The values of Pc 
and Ec are also indicated. 



nificantly smaller than Mtot indicating the disintegration 
of the shell into pieces. The value of the critical pressure 
Pc needed to achieve fragmentation and the functional 
form of the curve above Pc were determined by plotting 
{M,nax/Mtot) as a function of the difference \Pa — Pc\ 
varying Pc until a straight line is obtained in a double log- 
arithmic plot. The power law dependence on the distance 
from the critical point {Mmax/Mtot) ^ \Po ~ " for 
Po > Pc is evidenced by Fig.EJi), where a — 0.66 ± 0.02 
was determined. The finite jump of Mmax at Pc indi- 
cates the abrupt nature of the transition between the two 
regimes. For impact loading (Mmax/Mtot) proved to be 
a continuous function of Eo, however, it also shows the 
existence of two regimes with a transition at a critical en- 
ergy Ec- In Fig. Eh) we show (Mmax/Mtot) as a function 
of the distance from the critical point \Eo — Ec\, where Ec 
was determined numerically in the same way as Pc- As 
opposed to the pressure loading, (Mmax/Mtot) exhibits a 
power law behavior on both sides of the critical point but 
with different exponents (Mmax/Mtot) ^ \Eo ~ Ec\^ for 
Eo < Ec and (Mmax/Mtot) ~ \Eo - Eel"'' for Eo > Ec, 
where the exponents were obtained as /? = 0.5 ±0.02 and 
a = 0.66 ± 0.02. Note that the value of a coincides with 
the corresponding exponent of the pressure loading. 

We also evaluated the mass- weighted average fragment 
mass M as the average value of the ratio of the second 
M2 and first Mi moments of the fragment mass distri- 
bution subtracting the mass of the largest fragment from 
the sum . The behavior of M shows again clearly the 
existence of two regimes of the break-up process with a 
transition at the critical point Pc and Ec- Under pressure 
loading due to the abrupt disintegration M can only be 
evaluated above the critical point Pc, while for the im- 
pact case M has a maximum at the critical energy Ec 
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typical for continuous phase transitions. Fig.[3») demon- 
strates that in both foading cases M has a power law 
dependence on the distance from the critical point, i.e. 
M \Eo- -Bcl"^, and M ~ |Po - P^"^ with the same 
value of the exponent 7 — 0.8 ± 0.02. 

The most important characteristic quantity of our sys- 
tem which can also be compared to the experimental re- 
sults is the mass distribution of fragments F{m). For 
pressure loading F[m) can only be evaluated above Pc. 
Under impact loading for Eo < Ec we found that F{m) 
has a pronounced peak at large fragments indicating the 
presence of large damaged pieces. Approaching the criti- 
cal point Ec the peak gradually disappears and the distri- 
bution asymptotically becomes a power law at Ec. Above 
the critical point the power law remains for small frag- 
ments followed by an exponential cut-off for the large 
ones. For the purpose of comparison, a mass distribution 
F(m) obtained at an impact energy close to the critical 
point Ec, and distributions at two different pressure val- 
ues Po of the ratio 1.6 are plotted in Fig.|5]along with the 
experimental results. For impact an excellent agreement 
of the experimental and theoretical results is evidenced. 
For pressure loading, the functional form of F{m) has 
a nice qualitative agreement with the experimental find- 
ings on the explosion of eggs, furthermore, distributions 
at the same ratio of pressure values obtained by simu- 
lations and experiments proved to show the same ten- 
dency of evolution, see Fig. |21 Figs. Et, d) demonstrate 
that by rescaling the mass distributions above the criti- 
cal point by plotting F{m) ■ M as a function oim/ M an 
excellent data collapse is obtained with 5 — ± 0.03. 
The data collapse implies the validity of the scaling form 
F{m) ~ m""^ • f{m/M), typical for critical phenomena. 
The cut-off function / has a simple exponential form 
exp{—m/M) for impact loading, and a more complex one 
containing also an exponential component for the pres- 
sure case. The rescaled plots make possible an accurate 
determination of the exponent r, where r = 1.35 ± 0.03 
and T = 1.55±0.03 were obtained for impact and pressure 
loading, respectively. Hence, a good quantitative agree- 
ment of the theoretical and experimental values of the 
exponent r is evidenced for the impact loading of shells, 
however, for the case of pressure loading the numerically 
obtained exponent turned out to be higher than in the 
case of exploded eggs. 

In the fragmentation of bulk systems under appropri- 
ate conditions a so-called detachment effect is observed 
when a surface layer breaks off from the bulk and under- 
goes a separated fragmentation process jfi,, .£] . This effect 
also shows up in the fragment mass distributions in the 
form of a power law regime of small fragments of an expo- 
nent smaller than for the large ones. Our results on shell 
fragmentation can also provide a possible explanation of 
thisjtind of composite power laws of bulk fragmentation 
and might even be relevant for the study of mass 




I E^-Ejl (arb, units 



10 ' _10" 



FIG. 5: Normalized mass of the largest fragment a) and av- 
erage fragment mass b) as a function of the distance from the 
critical point. Rescaled plots of the fragment mass distribu- 
tions c) for explosion, and d) for impact. 



redistribution following super nova explosions in the uni- 
verse. 
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